Mobile tt— kinks and half-integer zero-field-like steps in highly discrete alternating 

— 7r Josephson junction arrays 
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The dynamics of a one-dimensional, highly discrete, linear array of alternating 0— and tt — Joseph- 
son junctions is studied numerically, under constant bias current at zero magnetic field. The calcu- 
lated current - voltage characteristics exhibit half-integer and integer zero-field-like steps for even 
and odd total number of junctions, respectively. Inspection of the instantaneous phases reveals that, 
in the former case, single tt— kink excitations (discrete semi-fluxons) are supported, whose propaga- 
tion in the array gives rise to the 1/2— step, while in the latter case, a pair of it— kink - 7r— antikink 
appears, whose propagation gives rise to the 1— step. When additional 27r— kinks are inserted in 
the array, they are subjected to fractionalization, transforming themselves into two closely spaced 
7T— kinks. As they propagate in the array along with the single tt— kink or the n— kink - n — antikink 
pair, they give rise to higher half-integer or integer zero-field-like steps, respectively. 

PACS numbers: 74.50.+r, 74.81. Fa, 75.10.Pq 



I. INTRODUCTION 

Conventional (0—) Josephson junctions (JJs) in the 
ground state, have zero 4>, the gauge invariant phase dif- 
ference of the superconducting order parameters on each 
side of the barrier. However, in some cases a J J may 
have a ground state with </> equal to tt (tt— J J). The pos- 
sibility of a 7T— JJ was first pointed out theoretically in 
the tunneling through an oxide layer with magnetic im- 
purities in SIS junctions, because of negative Josephson 
coupling [l[. Negative coupling also arises in JJs with 
ferromagnetic intermediate layers @, HJ . In that case the 
state of the J J, i.e., a 0— state or a tt— state depends sen- 
sitively on the sample design and the temperature 
Another realization of 7r— JJs became possible due to the 
d— wave symmetry of the order parameter in high— T c 
superconductors. Thus, tt— JJs were demonstrated using 
twist-tilt, and 45° tilt YBCO grain boundary JJs 
Very recently, a new technique of tailoring the barrier 
of superconductor- insulator- ferromagnet-superconductor 
(low T c ) junctions can be used for accurate control of the 
critical current densities of tt— JJs (as well as 0— JJs and 
— 7T JJs) @, [!]• All these junctions are characterized 
by an intrinsic phase shift of tt in the current-phase re- 
lation or, in other words, an effectivelly negative critical 
current. 

One-dimensional (ID) discrete arrays of parallel biased 
0— JJs for which both dimensions are smaller than the 
Josephson length Xj (i.e., small JJs), have been studied 
extensively in recent years 0, [13]. Such systems rep- 
resent an experimental realization of the spatially dis- 
crete sine-Gordon (SG) lattice, whose dynamic behaviour 
has been modelled with the discrete sine-Gordon (DSG) 



equation. That equation is formally equivalent to the 
Frenkel-Kontorova model, which describes the motion of 
a chain of interacting particles subjected to an exter- 
nal on-site sinusoidal potential [Tlj |. Experiments and 
numerical simulations have shown that, even with large 
discreteness, the dynamics of a localized kink (fluxon) 
in the SG lattice exhibit some features of solitonic na- 
ture, close to the properties of the continuous SG solitons 
j Til . However, discreteness also introduces a number of 
aspects in fluxon dynamics which have no counterparts 
in the continuum system. The experimental investiga- 
tion of discrete JJ arrays is a relevant issue in super- 
conducting electronics because they are the basis of the 
so called phase-mode logic [l3[ and of the rapid single- 
flux-quantum (RSFQ) circuits [3], for the construction 
of superconducting quantum interference niters (SQIFs) 
[l "J , as well as for developing high frequency mixers [l6| , 
or variable inductors and tunable filters [17] . 

Naturally, the next step in this field was to consider 
JJ arrays comprising tt— JJs. The usage of tt— JJs as 
passive tt— shifters in RSFQ circuits was recently sug- 
gested in reference [T^], and realized in pjj]. More- 
over, recent studies (both theoretical and experimental) 
of JJ arrays with tt— JJs have shown some novel fea- 
tures arisin g ou t of the interplay between 0— and tt— JJs 
[H EH, II M EE H, HI m E! H • F «r instance, 
when a 2tt kink is introduced in a hypothetical parallel 
array with alternating 0— and tt— JJs [31], it will break 
up into two separate tt— kinks. This effect is reffered to as 
"fractionalization of a flux quantum" and, surprisingly, 
it also appears in long continuous JJs with alternating 
and tt facets [32j |. In the limit of very short, randomly 
alternating and tt facets, that problem has been stud- 
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FIG. 1: Schematic view of a linear array of alternating 0— 
and n— JJs in parallel. Solid lines show the superconducting 
electrodes, while the crosses and the crosses- in- box show the 
0— and 7T— JJs, respectively. 



ied in references |33l l34j. Given that n— JJs can now be 
formed in a variety of ways, a parallel array of alternating 
0— and 7T— JJs could, in principle, be made in the lab. It 
is worth noting that parallel JJ arrays consisting of both 
0— and 7T— JJs were also studied as a model of long, 
faceted, high— T c bicrystal JJs with high misorientation 
angle of the bicrystalline film (~ 45°) Si, S3, SI, III . 

In the present work, we investigate numerically the 
fluxon dynamics in a ID, highly discrete array of alter- 
nating 0— and 7T— JJs. We are mainly concerned with 
linear arrays, with open ended boundary conditions, and 
different (even or odd) total number of JJs. In the next 
section, we describe the array model and typical ground 
states of such arrays, while in section 3 we present numer- 
ically obtained, typical current - voltage characteristics 
where zero field-like (ZF-like) steps appear. We finish in 
section 4 with the conclusions. 



II. MODEL EQUATIONS AND GROUND STATE 

Consider a ID array of alternatingly 0— and n— JJs, 
which are connected in parallel via superconducting 
leads, as shown schematically in figure 1. The Hamil- 
tonian function of this system in the simplest approxi- 
mation, where all mutual inductances between different 
cells are neglected, is given by 

n 

e"' { 2 k ^ n ~ ^ n ~ 1 ) 2 + t 1 ~ Icn cos ^") -7<M ,(1) 

where <j> n is the phase difference over the n— th JJ, I cn 
is the critical current of the n— th JJ, a is the dissipa- 
tion coefficient due to quasiparticles crossing the barriers, 
k = 1 /a 2 is the coupling strength between nearest neigh- 
bouring JJs, with a being the discreteness parameter, 7 
is the applied constant (dc) bias current assumed to be 
injected and extracted uniformly (see figure 1), and 



Pr, 



+at 



dt 



(2) 



is the canonical variable conjugate to <f> n . For a = 
that variable reduces to the (normalized) voltage - phase 



relation, v n — -ffi, for the n— th JJ. Since H depends 
explicitly on time t, it is obvious that it is not a constant 
of the motion. However, it is interesting that the system 
in consideration belongs in a certain class of systems hav- 
ing time-dependent Hamiltonians, which may be useful 
when studying stability issues of specific solutions. For 
a = and 7 = we get the simpler Hamiltonian 

+ ^ {(1 - I cn cos<j) n )} , (3) 



which has been used (with I cn = const.) to describe a 
harmonic chain of particles subjected to an external si- 
nusoidal potential In this context, the variable cj> n 
represents the displacement of a particle from its equilib- 
rium position, and the three terms in the earlier Hamilto- 
nian H s are interpreted as the kinetic energy of the par- 
ticles, the harmonic (clastic) interaction of the nearest 
neighbouring particles in the chain, and the interaction 
of the chain with an external on-site sinusoidal potential, 
respectively. If the ratio k of the coefficients of the elas- 
tic coupling energy to the on-site potential energy is less 
than unity, then the system is considered to be highly 
discrete. In respect to kink propagation, the system is 
highly discrete if the kink width is of the order of lattice 
spacing a. 

The n— dependent critical current I cn is given in com- 
pact form by 



= / c [£ + cos((n + l)7r)], 



where 



IcO ~ Ic 



(4) 



(5) 



with I c o (Ictt) being the critical current of the 0— (jr— ) 
JJs. From equations ^ and ([5]) we have that I cn = I c q 
and — I C7r for odd and even n, respectively. Note that I cn 
may differ both in sign and magnitude for the 0— and 
7T— JJs. From the Hamiltonian function given in equation 
((!]) we get the equations 



dt 2 

k{4>n+l 



d<f> n 

a—— 

dt 

24>n 4 



4>n- 



smt 



7. 



(6) 



For a finite array with N total number of JJs, n = 
1,2,...,N. The earlier equations may also be obtained by 
applying Kirchhoff laws in the equivalent electrical circuit 
[36(, where each single J J is described by the resistively 
and capacitively shunted junction (RCSJ) model. In this 
model the JJ consists of a capacitive branch, a resistive 
branch, and a superconducting branch, all connected in 
parallel. 

Equations ^ are written in the standard normaliza- 
tion form, i.e. a = D/Xj is normalized to Xj with D being 
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FIG. 2: The normalized flux 27r$ n /$o vs plaquette number 
n, for (a) a periodic array (iV = 100) ; (b) an open-ended array 
with even N (N = 100); (c) an open-ended array with odd N 
(N = 101). Ico = hn, and k = 0.8. 
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FIG. 3: Current-voltage characteristics for an array with 
N even (N = 100), I c o = Ic-w, ct — 0.2 and k as shown in 
the figure. The arrows along the dotted lines indicate the 
direction of the jumps. Inset: half-integer ZF-like steps in 
the lower-voltage part of the curves shown in main figure. 
The vertical arrows indicate the position of v m for ZFS1/2 
and ZFS3/2 for k = 0.8 (see text). 



the spatial separation of neighbouring JJs, the tempo- 
ral variable t to the inverse plasma frequency u>~ , and 
7 = Ic/IcQ is normalized to I c q. Then, I c q = 1 in equa- 
tions ©, and I C7r actually represents the ratio I C7r /I C Q. 
The voltage V n across the n— th JJ is normalized so that 
v n = V n /Vp where V p = <J>ow p /(27r) , and $o = h/(2e) is 
the flux quantum (with h and e being the Planck's con- 
stant and the electron charge, respectively). Equation 
is analogous to that deduced in reference [13] in the 
context of analysis of zig-zag arrays [38| . Below we deal 
with arrays both with periodic and open ends, i.e., annu- 
lar and linear finite arrays, respectively. Then, in both 
cases, the governing equations ^ remain the same for 
n = 2,,..,N — 1. The same equations are also valid at 
the end points n = 1 and n = N, if we define artificial 
phases </>o and </>jv+i such that 

0o (t) = <Pn (t) - 2irM f <p N+1 (t) = 0r (*) + 2%M > , (7) 
for periodic boundary conditions (e.g., [3(|[3!|), and 

Mt) = Mt) &v+i(t) = M*)> (8) 

for open-ended boundary conditions. The integer Mf in 
equation |[7J) is the number of 2tt fluxons trapped in the 
annular array. 

Equations ^ are integrated with a standard fourth 
order Runge - Kutta algorithm with fixed time-stepping 
(typically 0.01). In figure 2 we plot the calculated nor- 
malized flux $„/$o = {4>n — 0n+i)/(27r) in each elemen- 
tary cell (plaquette) as a function of the junction number 
n, for three different arrays in their ground state: a pe- 
riodic array (figure 2(a)), and open-ended arrays with 
even and odd N (figures 2(b) and 2(c), respectively). 
In the former case the flux clearly alternates between 
two values of opposite sign and the same magnitude, the 
latter depending on the coupling constant k, indicating 



that the flux lattice is actually periodic with a period of 
two elementary cells. We can thus define a new extented 
unit cell, with respect to the flux periodicity, which con- 
sists of two elementary cells. Then, we can evaluate the 
(normalized) average flux in each of the extended cells, 
V = ($2m-i+*2m)/(2$ ), where to = 1, N/2 num- 
bers the extended cell composed of the (2m — l)th and 
2mth elementary cells. Clearly, $ pm vs to is a constant 
which equals to zero. The same remarks also hold for the 
open-ended arrays, both with even and odd N, relatively 
far from the end JJs. 



III. CURRENT - VOLTAGE 
CHARACTERISTICS AND ZERO-FIELD-LIKE 
STEPS 

For obtaining the current-voltage (7 — v) character- 
istics for an open-ended array, we initialize the system 
with <f> n = and tt for n odd and even, respectively, and 
<p n = for any n at t — and 7 = 0. (The overdot de- 
notes derivation with respect to the temporal variable.) 
Then equations © are integrated forward in time up 
to t\ — 1000, by which time the system has reached a 
steady state. We then calculate the time averaged volt- 
age v n across each JJ by averaging over an additional 
t%— 1\ dimensionless time- units, i.e., v n = t ^_ t J t * 2 4>n dt, 
where £2 = 2000. The voltage of the array v is obtained 
by averaging over all JJs, i.e., v — J2n=i V n/N. The 
dc bias 7 is then increased by £7 and the calculation is 
repeated, assuming as initial condition the steady state 
achieved in the previous calculation. We continue this 
procedure increasing 7 by £7 until 7 = 1.1, then de- 
crease 7 with same steps back to zero, obtaining the 
7 — v characteristics shown in figures 3 and 4 for even 
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FIG. 4: Current-voltage characteristics for an array with 
A^ odd (N — 101), IcO = Ictt, ot — 0.2 and k as shown in 
the figure. The arrows along the dotted lines indicate the 
direction of the jumps. Inset: integer ZF-like steps in the 
lower-voltage part of the curves shown in main figure. The 
vertical arrows indicate the position of v m for ZFS1 and ZFS2 
for k — 0.8 (see text). 
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FIG. 5: An instantaneous profile of (a) $n/$o vs n, and (b) 
$ pm vs. extended cell number m, for J c o = I C tt, N — 100, 
7 = 0.2, k = 0.8, and 7 = 0.13 on the ZFSl/2-like branch in 
the inset of figure 3. 



2V(=100) and odd iV(=101), respectively (£7 = 0.005 for 
the curves in the main figures and 0.001 for the curves 
in the insets). Large hysteresis is observed in the high- 
voltage part (v ~ 0.5) of all characteristics in both fig- 
ures, even though the damping coefficient is rather high 
(a = 0.2). Small hysteretic effects are also apparent in 
the low- voltage parts (not shown for clarity). Appar- 
ently, the high-voltage characteristics in figures 3 and 4 
for same k are practically the same. The solutions on the 
high vo ltag e branches have been discussed in detail in ref- 
erence [3 II ] , and it seems that they are not significantly 
affected by the choice of boundary conditions. Here, we 
focus on the the low- voltage parts of those characteristics, 
which exhibit significant quantitative differences for odd 
and even N. That can be seen by comparing the insets 
of figures 3 and 4, which display a magnified portion of 
the corresponding low-voltage characteristics seen in the 
main figures. There we can see the appearence of current 
steps, similar with those observed in 0— JJ arrays (see for 
example references [IH, [HI). Those steps are the discrete 
counterparts of the zero-field steps (ZFSs) appearing in 
the 7 — u characteristics of continuous, long 0— JJ, due 
to fluxon resonant motion. In this case the propagating 
fluxons are 2ir— kinks and their multitudes, while single 
7T— kinks are unstable. In a continuous, long 0— JJ, with 
equivalent normalized length L ~ Na, the normalized 
voltage at maximum fluxon velocity v m of the first ZFS 
would be v m = 2n/L ~ 0.056 for k = 0.8. Thus, the first 
half-integer ZFS (ZFS1/2) would have v m = 0.028, while 
the second integer ZFS (ZFS2) and second half-integer 
ZFS (ZFS3/2) would have v m = 0.112 and 0.084, repec- 
tively. The fluxon width decreases (i.e., the fluxon size 
contracts) with increasing velocity. Thus, in the discrete 
system, a fluxon cannot reach the maximum velocity of 
its continuous counterpart, since its contracted size can- 



not be smaller than a. In each of the insets in figures 

3 and 4 there are two ZFS-like branches with k = 0.8, 
with maximum v (defined at the top of the steps) equal 
to 0.024, 0.072, and 0.048, 0.093, respectively. Thus, the 
steps in the inset of figure 3 can be identified with the 
discrete counterparts of half-integer ZFSs (namely with 
ZFS1/2 and ZFS3/2), while those in the inset of figure 

4 with the discrete counterparts of integer ZFSs (ZFS1 
and ZFS2). 

Inspection of the time evolution of the phases (not 
shown) reveals that the ZFSl/2-like branch for even 
Af(=100) is due to a ir— kink propagating in the array, 
on a staggered phase background. In figure 5 we show 
the instantaneous profile of the normalized flux as a func- 
tion of the plaquette number n, and $ pm as a function of 
the number of the extended cell m, for 7 = 0.13. The an- 
tiferromagnetic background in the flux (figure 5(a)) ob- 
scures the picture, but its proper elimination with the 
doubling of the unit cell along with the definition of $ pm 
makes things more clear (figure 5(b)). Then, in figure 
5b we see clearly the discrete n— kink in a form which is 
the discrete counterpart of the spatial phase derivative in 
continuous systems. Also apparent they are small oscil- 
lations behind the it— kink, a phenomenon known from 
0— JJ arrays [H, H^|, where moving 27r— kinks excite lin- 
ear oscillating modes of the array (Cherenkov radiation) . 
Those oscillations, however, get strongly damped and 
die out fast away from the kink. In figure 6 the same 
quantities for an array with the same parameters but 
odd AT(=101) are plotted for 7 = 0.13 on the ZFSl-like 
branch. Here we always have a tt— kink and a 7r— antikink 
pair. Cherenkov radiation is also apparent here (figure 
6(b)) behind both the propagating ir— kinks. Inserting 
one (or more) 2ir— kink(s) in both the even and odd 
N arrays, we get higher order ZFS-like branches. The 
ZFS3/2- and ZFS2-like branches for even and odd N, 
respectively, are shown in the insets of figure 3 and 4, 
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FIG. 6: A snapshot of (a) $„/$o vs n, and (b) $ pm vs. 
extended cell number m, for I c o = I C7r , iV = 101, 7 = 0.2, 
k — 0.8, and 7 = 0.13 on the ZFS1— like branch in the inset 
of figure 4. 



FIG. 7: A snapshot of (a) $„/<I>o vs n, and (b) <3> pm vs. 
extended cell number m, for I c o = /ctt, A = 100, 7 = 0.2, 
= 0.8, and 7 = 0.13 on the ZFS3/2— like branch in the inset 
of figure 3. 



respectively (only for k — 0.8). Those branches were ob- 
tained by inserting one 2n— kink in the arrays. Again, 
inspection of the phases reveals that the 27r— kink is frac- 
tionalized into two closely spaced but clearly separated 
7T— kinks which, however, are not identical. The snapshot 
in figure 7 shows the normalized flux as a function of the 
plaquette number n, and $ pm as a function of the number 
of the extended cell m, for 7 = 0.13 on the ZFS3/2-like 
branch in the inset of figure 3 (N = 100, even). The 
fractionalized 2n— kink at the right of figure 7(b), now 
coexists with the single 7r— kink in the array. The corre- 
sponding figure (figure 8) for the ZFS2-like branch in the 
inset of figure 4 (N = 101, odd), shows the fractionalized 
2ir— kink along with a tt— kink - n— antikink pair. 

The direction of motion of the fractionalized 27r— kinks 
depends on their polarity when they are inserted in the 
array. These kinks, either the single one or the ones 
obtained by fractionalization, behave in many aspects 
like the usual 2ir— kinks. They get reflected by the ends 
changing their polarity, and interact completelly elasti- 
cally, passing the one through the other without any ap- 
parent shape change. Although not shown in figures 3 
and 4, the depinning current 7 C related to the Pcicrls- 
Nabarro potential, which the kink should overcome in 
order to start moving, is practically the same for N even 
and odd. It increases with increasing discreteness (de- 
creasing k), having the values j c = 0.123, 0.099 and 0.066 
for k = 0.8, 1.0 and 1.5, respectively. The fractionalized 
27T— kinks shown for a specific illustrative case in the fig- 
ures 7(b) and 8(b), are closely spaced and they move 
together in the same direction without changing their 
distance of separation. Thus, they essentially form a sin- 
gle bound state which resembles the bunched fluxon state 
observed in 0— J J arrays [4(j. Such a state results from 
strong interactions between kinks, which is mediated by 
their Cherenkov "tails". As a result, the Cherenkov ra- 
diation behind the bunched state is strongly suppressed, 



which is exactly what we observe in our figures 7(b) and 
8(b). 
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FIG. 8: A snapshot of (a) $ n /$o vs n, and (b) $ pm VS. 
extended cell number m, for I c o = I C tt, N = 101, 7 = 0.2, 
k = 0.8, and 7 = 0.13 on the ZFS2— like branch in the inset 
of figure 4. 

In annular — 7r JJ arrays with both even and odd N, 
which satisfy the boundary conditions ([7]), the effect of 
fractionalization of 27r— fluxons also appears [3l|. In that 
case, 27T— fluxons can be inserted in the array either by 
setting Mf to a nonzero integer value or by an appro- 
priate initial condition. For both even and odd and a 
nonzero dc bias 7 (e.g., 7 = 0.13), there is a number of 
7T— kinks and tt— antikinks moving along the array such 
that the total flux is always an integer as required by 
the boundary conditions. They appear either as ir— kink 
- 7T— antikink pairs moving in opposite directions, or as 
two bunched 7r— kinks looking very much like those in fig- 
ures 7 and 8. However, we did not find single 7r— kinks 



6 



in this case, so that half-integer ZF-like branches do not 
appear. Thus, it seems that half- integer ZF-like branches 
appear only for arrays with even N satisfying the bound- 
ary conditions l[8"|). i.e., for linear arrays. This is probably 
due to the amount of the net total flux in a — tt J J array 
in its ground state. For arrays satisfying periodic bound- 
ary conditions without inserted flux, with either odd or 
even N, the ground state is always arranged such that 
the total flux is equal to zero. Apparently, this is not al- 
ways the case for open-ended boundary conditions. For 
odd N there is an even number of elementary cells in 
the array, so that the contributions from all cells exactly 
cancel, resulting in zero total flux. However, for even N 
there is an odd number of elementary cells in the array, 
and the total flux attains a nonzero value (approximatelly 
equal to 0.274 in figure 2(b)). With increasing bias cur- 
rent above the Peirls-Nabarro barrier, that ground state 
becomes unstable against the formation of a tt— kink. Dy- 
namic simulations for such — tt JJ arrays with periodic 
boundary conditions have been performed in reference 
[oil ], where it was shown that 2n— kinks fractionalize into 
a pair of 7r— kinks in certain regions of k — I cp space. 
However, ZF-like steps on the current - voltage charac- 
teristics are not discussed. Although a — tt JJ array 
is still hypothetical, it can be possibly realized with the 
existing technology. Moreover, it has been argued [13] 
that bicrystal grain boundary JJs with 45° misorienta- 
tions angle, form intrinsically an array of alternatingly 
and tt contacts. In this context, a discrete array of alter- 
nating 0— and 7T— JJs has been used as a model to study 
the dynamics as well as the magnetic field dependence of 
such JJs [13, Hi. 



allel connected alternating 0— and tt— JJs, which exhibit 
current steps in the low- voltage part. For even N, the 
maximum voltage of the steps (evaluated at the top of 
each step) is approximately two times smaller than the 
maximum voltage of the steps for odd N . Comparison 
with the continuous long J J, helps in identifying those 
steps as half-integer and integer ZFS-like steps for even 
and odd N, respectively. Inspection of the phases and 
the fluxes reveals that in both cases the system sup- 
ports 7T— kinks, either a single tt— kink or a tt— kink - 
tt— antikink pair for even and odd N, respectively, which 
can propagate in the array. Moreover, any 27r— kink in- 
serted in the array is fractionalized into two, generally 
different tt— kinks, which propagate together in the array 
along with the single tt— kink or the tt— kink - tt— antikink 
pair, for even and odd N, respectively, giving rise to 
higher order half-integer and integer ZFS-like branches. 
In the continuum, half-integer ZFSs have been observed 
experimentally in not very long — tt JJs due to the hop- 
ping of the semifluxons [41] ■ Similar ZFS-like branches 
can be found in rather wide parameter range. The case 
with k = 0.8, for which the system is already highly 
discrete, was chosen only for better illustration of the 
half-integer ZFSs. Those steps can also be found even 
for low damping (i.e., a = 0.02), as well as for higher 
damping (i.e., a — 0.42), which is more appropriate for 
high-temperature superconductor J J arrays [381 ] . Also, 
they may exist even for slightly different critical currents 
(i.e., for I C7T = 0.9 for even N). Dynamic states such as 
those described above could in principle be detected by 
low-temperature scanning electron microscopy (LTSEM) 
imaging [421 ]. 



IV. CONCLUSIONS 
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